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Abstract 

We show that the Lagrangian for interacting nonrelativistic particles can be coupled 
to an external gauge field and metric tensor in a way that exhibits a nonrelativistic 
version of general coordinate invariance. We explore the consequences of this invari- 
ance on the example of the degenerate Fermi gas at infinite scattering length, where 
conformal invariance also plays an important role. We find the most general effec- 
tive Lagrangian consistent with both general coordinate and conformal invariance to 
leading and next-to-leading orders in the momentum expansion. At the leading order 
the Lagrangian contains one phenomenological constant and reproduces the results of 
the Thomas-Fermi theory and superfluid hydrodynamics. At the next-to-leading order 
there are two additional constants. We express various physical quantities through 
these constants. 
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1 Motivation 



Symmetry is one of the most important ideas in modern physics. Frequently, the identifica- 
tion of the symmetries of a physical system is the key to understanding it. In this paper we 
identify a new type of symmetries of nonrelativistic physics, which are manifest when the 
system is put in curved space, yet less evident in flat space. 

The uncovering of these symmetries comes from our attempt to understand the properties 
of the so-called "unitary" Fermi gas, which is a gas of fermions which interact with each other 
through a zero-range two-body potential tuned to infinite scattering length a = oo. This 
system has been studied theoretically for some time [1, 2, 3], but only recently did it become 
realized in experiments with trapped atoms using a Feshbach resonance [4, 5, 6, 7, 8]. What 
is special about this system the absence of any intrinsic parameter, except for the number 
density. Unfortunately, this same feature means that there is no small coupling constant on 
which to build a reliable perturbation theory. It is therefore imperative that we find and 
exploit all the symmetries of this system. 

Despite the strong coupling, we are now certain that the ground state of the unitary Fermi 
gas is a superfluid. We also have, as our starting point, at least two theoretical frameworks 
which allow one to solve certain problems involving the unitary Fermi gas in a reliable fashion. 
These are the Thomas-Fermi theory and superfluid hydrodynamics. The two theories are 
very similar: both make no attempt to solve the problem at the microscopic level; they 
are both valid in the long-wavelength regime, when the characteristic length scale of the 
problem is much larger than all microscopic scales (including the healing length); and they 
need as an input at most one parameter (which is denoted as £ in Eq. (1) and characterizes 
the zero-temperature equation of state.) The total energy and the oscillation frequencies 
of unitary Fermi gas in a trap [9, 10] can be computed using these approaches. The input 
parameter £ is not fixed within these theories, but can either be measured experimentally or 
computed by Monte-Carlo simulations. 

However, until now it is not possible to systematically include corrections to either the 
Thomas-Fermi theory or superfluid hydrodynamics. At weak coupling there is no funda- 
mental obstacle to compute these corrections, but our tools do not work at strong coupling. 
These corrections can be considerable if we want to study the physics at length scales not 
much larger than the healing length. 

The main, utilitarian motivation of this paper is to develop a systematic way to go 
beyond the leading order in the long-wavelength expansion. This is done using an effective 
field theory (EFT) [11]. The central idea of the EFT philosophy is a power expansion 
over momentum (or inverse wavelength), instead of the usual perturbative expansion over 
powers of the coupling constant. This makes EFT particularly suitable for problems involving 
strongly coupled systems in the long- wavelength limit. Naturally we find that at leading order 
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the EFT reproduces the Thomas-Fermi theory and superfluid hydrodynamics. However, in 
principle, the EFT can be extended to arbitrary order in the expansion over momentum at 
the cost of introducing an increasing number of free parameters. The form of the effective 
Lagrangian is constrained only by the symmetries of the system under study. 

Somewhat unexpectedly, putting all possible constraints on the next-to-leading order 
(NLO) Lagrangian turns out to be a very nontrivial problem in the case of the unitary 
Fermi gas. In this paper, this task is accomplished by temporarily putting the system in 
curved space and exploiting the properties of general coordinate invariance and conformal 
invariance. Working in curved space necessitates the use of certain mathematical tools 
normally employed only in general relativity, but what we find at the end is quite nontrivial 
and worth the effort. In addition to one free parameter at the leading order, one needs just 
two parameters to characterize the NLO Lagrangian. These parameters cannot be computed 
within the EFT, but since the number of physical quantities that can be computed is greater 
than three, the EFT has real predictive power. 

We start the paper with a summary of some results that can be obtained using the 
EFT (Sec. 2). Then in Sec. 3 we show that the Lagrangian for interacting nonrelativistic 
particles can be coupled to background gauge fields and spatial metric in a way that respects 
gauge invariance and the general coordinate invariance. The general coordinate invariance 
is a gauge version of the translation and Galilean transformations. In Sec. 4 we study the 
conformal invariance of unitary Fermi gas. We then proceed, in Sec. 5 and 6, to find the 
most general effective Lagrangian, to leading and next-to-leading orders, that is consistent 
with these symmetries. Some applications of the EFT are considered in Sec. 7. We conclude 
with Sec. 8. 

2 Summary of results 

Here we give a brief summary of the result of this paper for the unitary Fermi gas. We 
do it here to give the reader a sense of the power and the limitations of the effective field 
theory derived in the paper. We work in the unit system where h — 1. We denote the 
number density as n and the fermion mass as m, from which we define the Fermi momentum 
k F = (3it 2 n) l l z , the Fermi energy e F = fcf/(2m), and the Fermi velocity v F = k F /m. We 
introduce the dimensionless parameter £, defined as the ratio of the energy density of the 
unitary Fermi gas to the energy density of a free Fermi gas at the same density, 

3 

e = £-e F n. (1) 
First we shall show that the dynamics of the gas is described, to leading order, by the 
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Lagrangian, 



£lo = com 



3/2 



H-V(t,x) -<p- 



2m 



5/2 



(2) 



where up is the phase of the condensate, \x is the chemical potential, V(t,x) is the trapping 
potential, and c is a dimensionless parameter. The leading-order Lagrangian captures the 
same information as the Thomas- Fermi theory and superfluid hydrodynamics. The coefficient 
Co is directly related to the parameter £, defined in Eq. (1), by 

2 5/2 



C 



(3) 



15tt 2 £ 3 / 2 ' 

and has to be determined microscopically. 

To the next-to-leading order in the low-momentum expansion, the effective Lagrangian 
has the form 

fVX) 2 



C = Com 3 ' 2 X 5 ' 2 + Cl m 1/2 - 



c 2 



'X \/i>) 

where we have introduced the shorthand notation 

X = fi-V(t,x)-<p 



(VV) 2 -9mV 2 A ]v / X. 



2m 



(4) 



(5) 



and ci and c 2 are two additional dimensionless parameters that appear at the NLO level. 

Many quantities can be computed using the effective Lagrangian. Let us consider infinite 
system first (V = 0). The dispersion relation for low-energy phonons is 



u(q) = c s q 



l-A/2£(c 1 + ^ 2 )!, 



k 



F J 



+ 0(q 5 lnq), 



The static density and transverse response functions are 



x(q) 



mkp 



l + 27rV 2 C ci 



9 \ r 



+ 0(g 4 lng), 



X T (q) = -9c n /^v F q 2 + 0(qnnq). 



(6) 



(7a) 



(7b) 



The predictive power of the EFT is illustrated by a relation between the nonlinear q 3 term 
in the phonon dispersion (6) and the static response functions, 

In addition, we also find the ground-state energy of a system of N particles in a harmonic 
trap to be 



E 



-co 



(3N) 4 ^-V2n 2 d Cl - 
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:°2 



Ujf + UJ 2 .+ 



(3iV) 2 / 3 + 0(iV 5 / 9 ), 



(9) 



where tUi, i = 1, 2, 3 are the trap frequencies, u = (cji^^) 1 ^ 3 - 
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3 Nonrelativistic general coordinate invariance 



3.1 Statement of general coordinate invariance 

Let us first consider a system of noninteracting particles (fermions or bosons). Instead of 
letting particles move freely in flat space, we will put the system in an external gauge field 

Ao = Ao(t,x), A i = A i (t,-x), i = 1,2,3. (10) 

and a curved three-dimensional manifold with the spatial line element 

ds 2 = gij (t, x) dx % dx^ . (11) 

As indicated in Eq. (11), the metric on the manifold is, in general, time-dependent, so the 
manifold can change shape with time. The spatial indices . . . will be raised and lowered 
by the metric and its inverse g x K The determinant of is denoted as g. The action of 
the system is 



S = J dtdxy/g C = J dtdxyfg , . ., . , 

(12) 

Here ip^ d t ip = ip^d t ip — d t ip^ip. Note that A , A iy and appear as external parameters of 
the Lagrangian; they have no dynamics of their own. The most common situation is when 
the system is put in a trapping potential V. This corresponds to A = V, Ai = 0, g^ = 5ij. 
Another case, not considered in this paper, is that of rotating systems where Ai are not zero. 
Although space is flat in all practical applications, keeping g^ arbitrary in intermediate steps 
allow one to derive nontrivial results, as we shall see. 

We have written the time-derivative term in Eq. (12) in the symmetric form dfip- 
In flat space this is equivalent to the asymmetric form iipW t ip; the two differ from each other 
by a full time derivative. This continues to be true if the metric is time-independent. For 
a general time-dependent metric, however, the two forms are not equivalent. As at the end 
we shall put the metric to flat, the choice between the two forms is a matter of convenience. 
We will soon find that the symmetric form in (12) is much more advantageous. 

The action (12) is invariant under gauge transformations, 

^_>^ = e ^, A -> A' = A - a, A % ^ A\ = A-fya, (13) 
where a = d t a. It is also invariant under a change of coordinate on the 3D manifold, 

V>(t, x) -> ^(t, x') = ip(t, x), A (t, x) -> A' (t, x') = Ao(t, x), (14) 

dx 1 dx^ dx l 
9ij(t,x) -> g V j>{t,-xT) = , , gij(t,x), A(i,x) -> A'(*, x = _^(t, x ). 
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The infinitesimal version of (13) and (14) is 



Sip = iaip - f k d k ip, 5A = -a - Cd k A , 

$9ij = -£ k dk9ij - 9ikdji k - g k jd£ k , 5A = -fya - £ k d k Ai - A k d£ k . 

We note that while the gauge parameter a in Eq. (15) can be a function of both time 
and space, so far we have considered the gauge parameter (which is the coordinate shift, 
x l — > x l + £ l ) to be time-independent. In fact, if one attempts to make £ l time-dependent, 
one finds that the action is not invariant, but acquires a change proportional to the time 
derivative of £ l , 

5S = "dtdx£ k ^d k ip. (16) 

However, the invariance of the action can be restored by a small modification of the trans- 
formation laws (15). The modified transformations read 

Sip = iaip - Z, k d k ip, (17a) 

SAo = -a- £ k d k A - A k £ k , (17b) 

5Ai = -dia - S k d k Ai - A k d£ k + mg lk £ k , (17c) 

tQij = -Z, k dk9ij - gikdjt k - g kj d£ k . (17d) 

It is straightforward to verify that terms proportional to £ k now cancel in SS. (Had we 
written the time-derivative term in the asymmetric form iip^d t ip this would not be true.) 
We shall call (17) the nonrelativistic general coordinate transformations, or simply general 
coordinate transformations. The action is said to be general coordinate invariant. 

Let us discuss the relationship between the transformations (17) with the more familiar 
Galilean transformations. In contrast to (17), the Galilean transformations are defined in 
flat space, ^ = Sij. They can be thought of as a special class of (17) with a particular choice 
of a and £ l : 

a = mwV, C = vH. (18) 
This leaves the flat metric unchanged while transforming other fields, 

ip(t, x) -> V>'(t, x) = e imv ^(t, x - vt), 

A (t, x) -> A' (i, x) = Ao(f , x - vt) - u fc A fc (t, x - vt), (19) 
Ai(t, x) -> AT(t, x) = x - vt). 

Note that in Eqs. (19) terms proportional to v 2 and higher orders of v have been dropped. 
Thus any theory with general coordinate symmetry is automatically Galilean invariant when 
restricted to flat space. The transformations (17) can be thought of as a gauge version of the 
translation and Galilean transformations. A similar "gauging" of nonrelativistic spatial 
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C = -tf d t ij - Aofy - ^ (d^ - iA^){d^ + iAjip) + q Q ^a - -g^d^a - ^ , (20) 



translations has been considered previously [12]. A subclass of (17) which does not change 
the flat metric has been considered in Ref. [13]. 

So far we have considered the noninteracting theory. It is easy to introduce interaction 
into the theory in a way that respects general coordinate invariance by introducing auxiliary 
fields. For example, the theory with the following Lagrangian, 

t •„ TW, , „ , _ _,. t _,._ 1 , _ 

is consistent with general coordinate invariance, as one can make a transform as 

5a = -i k d k a. (21) 

In flat space (20) describes a system of particles interacting though a Yukawa two-body 
potential, 

V(r) = -^ e - r/r °. (22) 
Anr 

We thus have an example of an interacting theory where coupling to an external gauge field 
and a metric can be introduced in a way consistent with general coordinate invariance. 

We have not demonstrated that general coordinate invariance can be achieved for any 
two-body potential. However, our foremost interest is in systems with universal behavior, 
where only the scattering length a is important, not the details of the two-body potential. 
For this purpose, the Lagrangian (20) is sufficient. In particular, infinite scattering length is 
achieved by fine-tuning g and r so that 

mqlr = 21.1... (23) 

when the first bound state appears at threshold. To achieve universality we also need to 
keep r smaller than any physical length scale. The universality of the fermionic system with 
infinite scattering length means that one can choose go and r arbitrarily (keeping r small), 
and the physics does not depends on our choice if Eq. (23) is satisfied. One can even make go 
and r time-dependent and nothing will depend on their values, if the fine-tuning condition 
is kept at all time. This remark will be important in our subsequent discussion. 

We have found the transformations (17) by trial and error. An alternative way to arrive 
to these transformations is by taking the nonrelativistic limit of the relativistic general 
coordinate invariance. The second way is instructive and will be useful for us when we try 
to write down an effective Lagrangian consistent with general coordinate invariance. We will 
describe it in the next Section. 

3.2 The relativistic root of nonrelativistic general coordinate invariance 

Consider a relativistic field theory of one free complex scalar field f in a external four- 
dimensional metric g^ u . We use the ( — h++) metric signature with x 11 = (ct, x) and ds 2 = 
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g ilv dx >x dx v '. We keep track of the speed of light c in order to send c — > oo in the nonrelativistic 
limit. The action, 

S = - Jd 4 x (g^d^d^ + m 2 c 2 ^) , (24) 
is invariant under the infinitesimal general coordinate transformation 



6* = -ed x y, 

$9^ = -£ x d x g^ - g\ u d^ x - g^\d u ^ x . 



(25a) 
(25b) 



In the nonrelativistic limit we take c — > oo. We assume that in this limit \1/ has a fast 
phase rotation in time solely due to its rest mass, which can be factored out to yield the 
nonrelativistic field ip, 

-imc 2 t 



1> 



= e 



\/2mc \/2mc 
The metric is taken to have the (00) component close to —1 and small (Oi) components, 



(26) 



9nv 



( 



\ 



2A 



A 

mc 



where A , A iy and gij are kept fixed in the limit c 



mc 

9ij J 

oo. The inverse metric is 



(27a) 
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(XV 



2A A l Ai 



m 2 c 2 



\ 



■— + 0(c- 3 ) 
mc 



A i \ 
■— + 0(c- 3 ) 
mc 



gV +0 (c- 2 ) J 



(27b) 



where A 1 = g tj Aj. Substituting Eqs. (26) and (27b) into the action (24), one reproduces 
(12) after discarding terms containing negative powers of c. Note that the goo and go% 
components of the 4D metric tensor now appear as the gauge field in the nonrelativistic 
action, in a manner reminiscent of the Kaluza-Klein mechanism [14]. 
Taking the coordinate shift ^ of the form 

/ a 



mc 



(28) 



where a is fixed in the limit c — > oo, Eqs. (25) reduce to Eqs. (17). Thus we can understand 
(17) as a nonrelativistic limit of the relativistic general coordinate invariance. 



3.3 General coordinate invariance and conservation laws 

We now show that the general coordinate invariance implies conservation of particle number 
and momentum. We also will find that in theories with general coordinate invariance the 
number current coincides (up to a factor of m) with the momentum density. 
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The connected part of the generating functional W is formally defined as a path integral 

e iW[A ,A t ,g t:i ] = / JJ^JJ^ e iS[ip^, A ,Ai,g i:i ] _ ^9) 



The correlation functions of operators coupled to the sources can be found by differentiating 
W with respect to its arguments. General coordinate invariance implies, to linear order in 
at and £ l , 

W[A + 5A , Ai + 8Ai, g tj + 5g tj ] = W[A ,Ai, 9ij ], (30) 

where 5A , 8A h and 5gij are defined as in Eqs. (17). Expanding the left-hand side of Eq. (30) 
to linear order in a and £\ one finds 

d t n + dif = 0, (31a) 
d t T ok + diT k = 0, (31b) 

where 

SW SW . SW 

n = -TT~ ' T ok = -mg ik —- + A k —- , 

oA dAi 6A 

$W . 5W T . Tr SW [ ' 

Eqs. (31) can be interpreted as conservation of particle number and momentum. Here n 
is the number density, j l is the number flux, T 0i is the momentum density and TJ, is the 
stress tensor. We also notice that in flat space and in the absence of the vector potential the 
momentum density is proportional to the number current, 

T 0i = mf when = 0, g tj = 5^. (33) 

This is the condition previously used by Greiter, Wilczek, and Witten [15] to find the effective 
field theory of superfluids. Any theory which is general coordinate invariant automatically 
satisfies (33). 



4 Scale and conformal symmetries 

The Lagrangian of noninteracting particles (12) exhibits scale invariance, which is an overall 
rescaling of time and fields, 

= \-H, 

x) - x) = A 3 /V(*, x), Ao(t) ^ A 'o(t') = \A (t), (34) 

9a(t) - qW) = Mt) - m>) = A t (t). 

It is useful to introduce the notion of scaling dimension. An operator O is said to have 
scaling dimension Ao (in shorthand notation, [O] = Ao) if under scale transformation in 
transforms as 

0(t, x) -> 0'(\-H, x) = X Ao O{t, x). (35) 
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Note that the scaling dimension of a product is the sum of scaling dimensions: [0]0 2 ] = 
[Oi] + [0 2 ]. Also, one can assign scaling dimensions to space and time derivatives, [<9j] = 
and [dt] = 1, so [dtO] = [didjO] = [O], while [d t O] = [O] + I. In this language, 

[t] = -l, [x i ]=0, ty] = l, [Ao] = l, [A] = 0, \9ij] = -l. (36) 

and [C\ — §, so S is invariant. 1 

In fact, the action is invariant under a larger set of transformations, which involves 
arbitrary reparametrization of time t — > t' — t'(t), simultaneous with rescaling of A and gij, 



(37) 



m - = m, Mt) - = (§) Mt), 

9ij(t) - gW) = (f ) 9ij(t), Mt) - A[(t') = Mt). 

The infinitesimal version of (37) is 

t -> f = t + 0, 

6^ = -(3i) - 5A = -PA - $Ao, (38) 

Sgij = -pgij + pgij, 5Ai = -fiM 

We shall call this transformation the conformal transformation. 2 The scale transforma- 
tion (36) is a particular case of conformal transformations with (3 being a linear function of 
time: (3 = bt. 

In general, if a field O transforms under conformal transformation as 

SO = -(30 - A $0, (39) 

then O is said to have conformal dimension A G , [0] c = A Q . Note that we still have [Oi0 2 ] c = 
[Oi] c + [02]o an( i [diO] c = [0] c . However, we cannot assign a conformal dimension to the 
time derivative d t . In fact, if O transforms as in Eq. (39), then in general O does not have 
a well-defined conformal dimension at all, 

SO = -(30 - (A + l)(30 - A (30. (40) 

The extra term proportional to (3 prevents us from assigning the conformal dimension Ao + 1 
to O, unless A Q = 0. (This term is absent for scale transformations where (3 is linear in t.) 

1 The zero scaling dimension of x and Ai in Eq. (36) may seem somewhat unnatural, but that is because 
we have assigned a scaling dimension of —1. 

2 In the literature [16, 17] the term "conformal transformation" was used for a particular a combination 
of (38) and the general coordinate transformation (17) with f3 = —Ct 2 , £ k = —Ctx k and a = -|mCx 2 , 
where C is a constant. This combination leaves the trivial background A = Ai = 0, = 5ij unchanged. 
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Therefore while an operator with a conformal dimension has a scaling dimension, the reverse 
is not true in general. 

Note that the Lagrangian (12) contains a term with a time derivative but still have a 
well-defined conformal dimension of |, since the j3 terms cancel out between the variations 
of ip^dfip and d t ip'ijj. Again we see that the symmetric from of the time-derivative term is 
preferable. 

The conformal symmetry is generally broken when one turns on an interaction. We now 
show that conformal invariance becomes exact again in the limit of infinite scattering length. 
For that we can use the model with Yukawa interaction (20). If we take a to transform the 
same way as ip, 

a(t) - a>(t>) = ^-J a(t), (41) 

then after substitution, we see that two terms in the action, q^ipa and — <r 2 /(2rl), change 
under the conformal transformation. The effect is that the constants q and r are replaced 
by new time-dependent constants 

However, we recall that the condition of infinite scattering length corresponds to mq^ro 
equal the critical value (23). We can see that this condition is not violated under the 
conformal transformation. Thus we can make use of the freedom in choosing q and r 
to transform it back to the original constant q and r , without changing the generating 
functional W. Therefore, the generating functional of a system tuned to infinite scattering 
length is conformally invariant. This property of unitary Fermi gas will be used in Sec. 6. 

Finally, we note that in contrast to the general coordinate transformations, the scale and 
conformal transformations described above are not descended from the relativistic counter- 
parts. 

5 The superfluid effective theory 
5.1 Preliminaries 

It is now well established that Fermi gas remain superfluid in the whole range of the pa- 
rameter kpa from the BEC to the BCS regimes. Superfluids provide the simplest nontrivial 
application of general coordinate invariance. At long distance the only physical degree of 
freedom is the phase of the condensate. We shall therefore construct an effective field theory 
of the this phase. 
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For a system consisting of bosons, the field is defined as the phase of the microscopic 
field, 

i> = He"**. (43) 
For fermions, we define 9 as half of the phase of the Cooper pair, 

{^) = \^)\e~ 2W . (44) 

Frequently, one describes a system at finite chemical potential /i by including the term 
/iip^ip into the Lagrangian. This approach is not suitable for our purposes since this term 
breaks general coordinate, scale and conformal symmetries. Instead we shall write the La- 
grangian without the [iip^ip term. The chemical potential enters the problem through the 
choice of the ground state. In particular, the superfluid ground state of a system at chemical 
potential ji corresponds to 9 = fit. We will frequently expand 

9 = fit-(p, (45) 

where fluctuations of ip around zero correspond to phonon excitations. The breaking of 
symmetries is now spontaneous, i.e., it comes from the choice of the ground state (45). 

Let us first discuss the theory in free flat space (A — Ai — 0, = 5^). The effective 
field theory is formulated as a Lagrangian L c r(9). This Lagrangian has a shift symmetry 
9 — > 9 + const, due to the global symmetry with respect to phase rotation of if). Due to 
this symmetry, 9 does not appear without derivative. So when we expand the Lagrangian 
in power series over derivatives and over 9, we encounter in each term at least as many 
derivatives as 9. Consider a few representative term in this expansion: 

C cS {9) ~ {89f + {89f + {8d9f + ■■■ (46) 

The most naive power counting scheme is to consider 9 ~ 1 and each derivative as 0{p). 
In this scheme, the first term in (46) is 0(p 2 ), while the last two terms are 0(p 4 ). However, 
the power counting scheme that is most useful is to count d t 9 and 8$ as 0(1), and additional 
time and space derivatives as 0{p). This is possible since, as we have seen, 9 always appears 
with at least one derivative. In this second power counting scheme, the first two terms in (46) 
are of order 1 while the last term is 0(p 2 ). A generic term in the which contains N[dt] time 
derivatives, N[di] spatial derivatives, and the N[9] power of the field will be counted as 

d N[d t ] g N[d t ] e N[e] ^ p N[d t ]+N[di]-N[e]_ ( 47 ) 

In the second power counting scheme, the leading order effective Lagrangian sums up all 
terms where the power of 9 is equal to the total number of space and time derivatives. Thus 
it is a function of first derivatives of 9, 

C LO = L(9,d t 9), (48) 
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which may contain, upon series expansion, terms with any power of dip. 

So far we have not discussed the background fields A , A^, and p^-. Since A and Ai 
appears with in covariant derivatives, it is natural to assume that all fields scale as 0(1) 
and their first derivatives as 0(p). In fact, in a trap where one has to treat the trapping 
potential A as a quantity of order 1, the second power counting scheme with 86 ~ 1 is the 
only consistent one. 

5.2 Symmetry requirements on the effective Lagrangian 

The effective theory should inherit the general coordinate invariance of the microscopic 
theory. This means that the effective action is invariant under the general coordinate trans- 
formations. Since 

$V9 = \V9 9 tJ 5g l3 = -fhy/g ~ y/gdk£ k , (49) 
we need the Lagrangian density to transform as 

5C = -£ k 8 k C, (50) 

so that the action is invariant 

5S = J dtdx {-i k C8 k ^ - y/gCd£ k - V9^ k dkC) = -J dtd^8 k {^i k C) = 0. (51) 

Equation (50) means that C has to transform as a scalar under general coordinate transfor- 
mations. 

We shall assume that the microscopic theory is invariant under time reversal, which 
means that the effective theory is invariant under 

t -> -t, 6^ -6. (52) 

This time reversal invariance has an important consequences for the momentum expansion 
of the EFT. This invariance implies that in each term, the difference between of the number 
of time derivatives and the number of 6 is even, 

N[8 t ] - N[6] = even. (53) 

Due to rotational symmetry and parity (which forbids e ufc ), the number of spatial derivatives 
have to be even, since the spatial indices have to contract pairwise. Therefore we have 

N[8t\ + N[8i] - N[6] = even, (54) 

which means that the terms in the effective Lagrangian are proportional to even powers of 
p. Therefore the next-to-leading order terms in the Lagrangian are 0(p 2 ), the next-to-next- 
to- leading terms are 0(p 4 ), etc. 
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5.3 Leading-order effective Lagrangian 

Let us now turn on and the curved metric. The leading-order effective Lagrangian, as 
we have discussed, should be a function of a finite number of variables which are 0(p°): 

C = C{b,d i 0,A Q ,A i ,9i3). (55) 

Gauge invariance and invariance with respect to three-dimensional general coordinate trans- 
formations (with time-independent £*) limit the Lagrangian to be a function of two variables, 

C = C (D t 9, ^DiODjB) , (56) 

where 

D t 9 = 9-A , D i = d i 6-A i . (57) 

Now we require the invariance of the effective theory with respect to general coordinate 
t r ansf or mat ions 

66 = - a - £ k d k e, 5A = -a- Z k d k A - A k ^ k , _ 

$9ij = -i k d k gij - 9ikdji k - g k jdit k , SAi = -dta - t h d k Ai - A k d£ k + mg ik £ k . 

Under (58) we have 

6(D t 6) = -i k d k D t 9 - i k D k 6, (59) 
8{g ij D i 9D j e) = —£ k d k (g ij Di9Dj9) - 2mi k D k 9. (60) 

This means that only one combination of these two variables transforms as a scalar, which 
is 

X = D t 9 - ^-D t 9Dj9, 5X = -^ k d k X, (61) 
and the leading-order effective Lagrangian has to be a function of this combination, 

C = P(X). (62) 

If we now set the metric to flat and external fields to zero, then we find the most general 
form of the Lagrangian for superfluids, 

C=p(»-®£), (63) 

which was previously found by Greiter, Wilczek, and Witten [15] using a different line of 
arguments. 
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Lagrangian density and thermodynamic pressure 

We now show that the function P in the leading-order Lagrangian is identical to the 
function that defines, at zero temperature, the pressure as a function of the chemical poten- 
tial, P = P(u). For that purpose, we notice that the particle number density can be found 
by differentiating the Lagrangian with respect to A Q , and is equal to 

n = P'(X). (64) 

The ground state at a chemical potential /i corresponds to 9 = fit, which means X = fi. This 
implies P is a function satisfying n(/i) = P'{fi) where n(/z) is the equilibrium number density 
at chemical potential /i. That means P is the thermodynamic pressure, up to a constant. 
The constant is not important in any situation and will be set to zero. 

Nonrelativistic superSuid as a limit of relativistic superRuid 

It is also possible to obtain the Lagrangian (63) by taking the nonrelativistic limit of the 
effective Lagrangian of relativistic superfluids. The field appearing in the effective Lagrangian 
is the phase 9 of the relativistic field \1/ (see Eq. (24), 

^ = \^\ e - ie . (65) 

From Eq. (26), the relation between 6 and the nonrelativistic phase 9 is 

6 = mch + 9. (66) 

To leading order, the effective Lagrangian is a function of <9 M 0. Due to Lorentz invariance, 
it has to have the form 

C = Co(T), T = -^d,Qd u e. (67) 



From Eqs. (66) and (27b) one finds 

+ 0( C - 2 ) = '-^—+mX + 0(c- 2 ). (68) 



^ m 2 c 2 

T = — — + m 



e - Ao- ^fre - AMe - Aj) 



lire 2 , _ 2 . 



So, in the nonrelativistic limit the Lagrangian is a function of X, as found in Eq. (62). 

It is shown in Ref. [18] that the function Cq{T) coincides with pressure at chemical 
potential equal to \ / 2c 2 T = mc 2 + X + 0(c~ 2 ). This is consistent with our previous result, 
since the difference between the relativistic and nonrelativistic normalizations of the chemical 
potential is the rest energy mc 2 . 

5.4 Relation to other approaches 

For completeness, we include here a discussion of the relation of the Lagrangian (63) and its 
derivation with other, more frequently used, approaches. We also derive some well-known 
results using the effective Lagrangian. 
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General coordinate invariance versus Galilean invariance 



Galilean invariance is sufficient to obtain the leading order Lagrangian. We recall that 
the Galilean transformations (18) are particular cases of general coordinate transformations. 
The Galilean transformations leave flat metric unchanged. Under these transformations, 

D t 9 -> D t 9 - ViDrf, Dfi -> Did - mvi, (69) 

from which we see that the only combination which is rotationally and Galilean invariant is 
D t 6-{2m)-\D i e) 2 . 

A natural question arises: does Galilean invariance in general guarantee general coordi- 
nate invariance? The answer is no. We shall see in Sec. 6 that at the next-to-leading order, 
Galilean invariance is not sufficient to constrain the form of the effective Lagrangian. 

Effective Lagrangian and superfluid hydrodynamics 

From Eqs. (63) and (45) we can derive the field equation for <p, 

d t n+-d i (nd i (p)=0, (70) 
m 

In Eq. (70) n is the equilibrium number density at chemical potential equal to X = 9 — 
(di9) 2 1 (2m) . This is the continuity equation of superfluid hydrodynamics, if we identify the 
superfluid velocity with the gradient of 9, 

v. = -™=^. (7!) 
m m 

In superfluid hydrodynamics one regards n and ip as independent variables. The second 
equation of superfluid hydrodynamics describes time evolution of ip comes from the definition 
of n, 

, . mv 2 

V = -Ufa) - ■ (72) 
Here /i(n) is the chemical potential at which the number density is equal to n. 

Phonon velocity and interaction 

We can determine the phonon spectrum by expanding the Lagrangian (63) to quadratic 
order in ip (recall that 9 = fit — ip). The Lagrangian is then 

£ = P M _^ + y>y_jL (M » (73) 

The first term is constant, and the second term is a full derivative and does not contribute 
to the field equation. The two last terms determine the phonon speed to be 
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In performing the last transformation we have use the thermodynamic relation dP = nd/j, 
and the definition of mass density p = mn. The speed of sound is given by the same formula 
as in normal fluids (however here the sound is not collisional.) 

Expanding the Lagrangian further, one can find the interaction vertices. The interaction 
between phonons, at leading order in momentum, is determined solely from the equation of 
state. This nontrivial result is actually not new. Its relativistic equivalence was derived in 
Ref. [18]. In the nonrelativistic setting, it was known in the Hamiltonian formalism, which 
we discuss next. 

Effective Hamiltonian 

From the effective Lagrangian (63) one can derive the equivalent Hamiltonian. The 
variable conjugate to the phase 9 is the number density, 

n^^ e = -.=P\X). (75) 
ov 

The equal-time commutation relation between n and 9 is 

[n(*,x),0(t,y)] = -i<5(x-y). (76) 

The Hamiltonian is then 

H = n9-P(X) =n^p^ + (nX - P(X)). (77) 

Note that the function nX — P(X) is the Legendre transform of P{X). Since P(X) is the 
pressure as function of chemical potential, its Legendre transform is energy as function of 
number density e(n). Therefore, 

H = n ^¥- + e(n). (78) 
2m 

The Hamiltonian is the sum of the kinetic energy pv^/2 and the internal energy not associated 
with the macroscopic flow, e(n). This Hamiltonian is well known. It was used, e.g., to find 
the interaction between phonons [19]. The Lagrangian — n9 — H is sometimes called Popov's 
hydrodynamic Lagrangian [20]. 

5.5 Effective Lagrangian for a unitary Fermi gas 

The equation of state of a unitary Fermi gas is known up to a constant. Due to the absence 
of any internal scale, by dimensional analysis the pressure has to be of the form 

P = c m 3 / 2 // 5 / 2 , (79) 
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where c is a dimensionless variable. The leading-order effective action is then, 

5/2 



S = J dtdxc m 3/2 (r> t e - ^-DiBDjO^ . 



(80) 



The fractional power in Eq. (80) is not problematic, since we always expand around the 
ground state 6 = fit. As expected, the Lagrangian (80) is invariant under infinitesimal 
conformal transformations 

59 = -/36, _ 6A = -(3A -$A , 

$9ij = -PfJij + P9ij: (5 'Ai = -13 Ai. 

Taking the Legendre transform of P([J>), we find the energy density as a function of the 
number density. This is proportional to the energy density of a free Fermi gas at the same 
density, with the coefficient of proportionality £, 

3/2n\ 2/3 n 3 (37r 2 n) 2 / 3 , . 

e=- — — = » — n - 82 

5 V5c y m 5 2m K J 

The parameter ^ is estimated to be around 0.4 from quantum Monte Carlo simulations [21, 
22] (but see also Ref. [23]). We can express c via £, 

2 5/2 

C ° = 15^ • (83) 

At large, but finite, scattering length a the effective Lagrangian contains an extra term 
proportional to a -1 , 

C = c m 3 / 2 X 5 / 2 + d — (84) 

a 

where d is some constant. The term proportional to d violates scale and conformal sym- 
metries explicitly. 



6 Superfluid EFT to NLO 

6.1 NLO terms versus loops 

The leading-order effective Lagrangian allows one to compute all quantities of interest, e.g., 
scatterings amplitudes or correlation functions, to leading order in momentum. The correc- 
tions to these leading-order results come from two sources. First, the Lagrangian contains 
terms subleading in derivatives; these terms contribute, at tree level, subleading contributions 
to the scattering amplitudes and correlation function. As explained above, these corrections 
start at 0(p 2 ) order. The second source of corrections are loop diagrams, constructed using 
the leading-order Lagrangian. 
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We shall demonstrate here that the loop contributions start at order 0(p 4 ), i.e., are 
comparable to those coming from the next-to-next-to-leading order terms in the effective 
Lagrangian. Therefore, if one is interested in the leading 0(1) and next-to-leading 0(p 2 ) 
results, one can work at tree level. 

To estimate the size of loop corrections, it is useful to work in the unit system where the 
sound speed is c s — 1, so energy and momentum has the same dimension. The dimension 
of the Lagrangian is 4 (momentum 4 ). We also rescale tp so that the quadratic terms in the 
Lagrangian is normalized canonically. The dimension of tp is then 1. The interaction terms 
will have coefficients which can be estimated from dimensional analysis, 



Here \i is the chemical potential and all dimensionless coefficients are, in general, of order 
one. 



As an example of a loop diagram, consider a one-loop correction to the phonon dispersion 
relation (Fig. 1). The graph is divergent and has to be regularized. After dimensional 
regularization it becomes a finite function of the external momentum p, with a possible 
logarithmic singularity. However we know that the graph contains two three-phonon vertices 
and hence is proportional to /i -4 . The self-energy has dimension 2. Therefore we conclude 
that 




(85) 




Figure 1: A one loop contribution to the phonon self-energy. 



Si-ioop jP 2 hip. 



(86) 



which is suppressed by 0(j9 4 lnp//x 4 ) compared to the tree-level result 0{p 2 ). 




Figure 2: A one loop contribution to the phonon scattering amplitude. 
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Consider another example, which is a one-loop correction to phonon-phonon scattering 
amplitude (Fig. 2). The scattering amplitude is 0(p 4 /fi 4 ) at tree level. The one- loop graph 
contains two four-phonon vertices and hence is proportional to /i~ 8 . By dimension, it has to 
behave as p s In p/fi 8 as a function of the external momenta p. The result is by a factor of 
0(p 4 lnp/ fi 4 ) smaller than the tree level. 

Generalizing to more complicated diagrams and processes, we see that we have to compute 
loops only when we want to know the answer to order p 4 compared to the leading-order 
result. If we want to know the answer only to the leading and next-to-leading order, tree- 
level calculation using the Lagrangian with subleading terms should suffice. 

6.2 Structures of NLO terms 

To find all terms allowed to the next-to-leading order, it is most convenient to work in the 
relativistic theory and then take the nonrelativistic limit. 

In the relativistic theory, the next-to-leading terms contain two derivatives more than 
the field 0. They have the form Q i F{X) ) i = 1 . . . 7, where F is an arbitrary function of X, 
and Qi are Lorentz scalars which are of order p 2 : 

Q 1 = (V^d u Q){V x d p Q)d^Qd' / Qd x Qd p Q = d^Td x Td^Qd x Q, 

Q 2 = (v^exne^we = -^r^ene, 

Q 3 = (v^exv^e^'we = d.Td^r, (87) 
Q 4 = (ne) 2 , q 5 = (VjAexv^e), q 6 = i^e&e, q 7 = r. 

Here V M is the covariant derivative, DO = V M <9 M 6, R^ u is the Ricci tensor, and R = R^ (see 
Ref. [24] and Appendix A). The terms proportional to the Ricci tensor are 0(p 2 ) since R^ v 
contains second derivatives of g^ u ~ 0(1). 

We now use the field equation to exclude redundant terms from (87). To leading or- 
der, where the Lagrangian is C = Cq(T), the field equation can be written as d^Td^Q = 
— (£'/£")□©; after adding 0{p 2 ) terms (87) to the Lagrangian, there will be higher correc- 
tions to this equation: 

^T^0 = -^D0 + O(p 3 ). (88) 
By using the field equation one can express Qi and Q 2 through Q 4 : 

Qi = (jp) Q± + 0( P 4 ), Q 2 = ^Q 4 + 0(p 4 ), (89) 
so they are not independent. In addition, one term from (87) can be written as a linear 
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combination of other terms, up to a full derivative, 

V At ^0V M ^0F(T) = (D0) 2 (f - ^j-F^j + d^Td"TF'(T) - R^d^Qd u QF(T) 

- V„ [(0"T + □G)F(T)] . (90) 

Therefore there are only four independent terms remaining from the original seven, and 

£nlo = d^Td^TFi (T) + (n&) 2 F 2 (T) + R^ed u QF 3 (T) + RF 4 (T), (91) 

where F X) F 2 , F 3 , and F 4 are arbitrary functions. 

Now all we have to do is to take the nonrelativistic limit of the terms appearing in (91). 
Consider the Fi term. Notice that 



d T = d {^-f + mlj = ™d t X = Oic- 1 ), 



(92) 



therefore 

d"Td^T = m 2 g ij d i Xd j X + 0{c~ 2 ). (93) 
For the F 2 term, we first recall that 

d"0 = g^d u e « (-mc, g ij Dj6), (94) 

therefore 

□0 = d ^ gg ^ dv Q) = dt y/g + —a i (g'>D j 0). (95) 

The coefficient in front of F 3 can be expanded as 

R llv & t Q&'Q = m 2 c 2 R 00 - 2mcR 0i D i 6 + R i jD i QD^Q (96) 

where all three terms survive in the limit c — > oo since i?oo = 0(c~ 2 ), i?oi = 0(c _1 ), and 
-Rjj = 0(1) [see Appendix A where the components of the Ricci tensor are computed for the 
nonrelativistic metric (27)]. 

The coefficient in front of _F 4 is the Ricci scalar, which reduces to the Ricci scalar of the 
3D space in the nonrelativistic limit after neglecting terms of order c~ 2 and higher: 

R = g^R^ w g ij Rij = R 3D (#,-)■ ( 97 ) 

Thus, we arrive to the most general NLO Lagrangian consistent with general coordinate 
invariance, 

£nlo = A + C 2 + C 3 + £ 4 , (98) 



21 



where 



C 1 =g ii d i Xd j Xf 1 {X), (99) 
C2 = - [-mdty/g + diiy/gg'WjetfMX), (100) 

f in 2 i 

C 3 = \ - — QgVgy + Fga) ~ m ^i E% + 7^ 

(101) 



+ 

C 4 = R m f 4 (X), (102) 

where E 1 , = d t Ai — diA , F i3 = c^A,- — <9jAj, and raising/lowering indices and covariant 
derivatives are done using the 3D metric g^. 

Now having made use of the general coordinate invariance, we dispense with the curved 
space and set the metric to flat, gij = Sij. The NLO Lagrangian becomes much simpler, 

Cnlo = diXdiXf^X) + {diDie) 2 f 2 (X) + (-mdiEi + -F tJ F l3 - d^D^ f 3 (X). (103) 

The fi term is not important in flat space. If one retains only A and sets A i = (which 
would preclude one from calculating the transverse response function, see Sec. 7), then the 
Lagrangian is even simpler 

£ NLO = diXdiXf^X) + (V 2 6) 2 f 2 (X) + mV 2 A f 3 (X). (104) 

General coordinate versus Galilean invariance 

We have mentioned above that, at the leading order, general coordinate invariance is not 
more powerful than Galilean invariance: both yield the Greiter-Wilczek-Witten Lagrangian. 
However, at the next to leading order, general coordinate invariance is more restrictive than 
Galilean invariance. In other worlds, one can write down Galilean invariant terms which are 
not general coordinate invariant. 

To see that, we write down the infinitesimal Galilean transforms (recall that we are in 
flat space), 

Dte^DtO-ViDiO, Ei^Ei-VkFu, 
D$ -> Dfi - mvi, F i:j -> F i3 . 

All three terms in Eq. (103) are Galilean invariant, as they should be since general coordinate 
invariance implies Galilean invariance. However, one can directly check that each of the 
following three expressions is Galilean invariant, 

FijFij, mdiE, + diFijDjO, m 2 E 2 + 2mE i F lk D k 6 + F i3 F lk D 3 QD k Q, (106) 
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but, in fact, the third expression does not enter the Lagrangian (103) at all, and the first 
two enter the Lagrangian only through one particular linear combination. Therefore if one 
relies only on Galilean invariance one would conclude incorrectly that the NLO Lagrangian 
involves five arbitrary functions, instead of three. When one only keeps the trapping potential 
A , setting = 0, general coordinate invariance allows one to exclude a term of the form 
(diA ) 2 f(X) from appearing in the Lagrangian (104). 

6.3 Scale invariance and conformal invariance 

Unitary Fermi gas does not possess an intrinsic scale, therefore the functions fi(X), i = 
l,---4 above should be powers of X. Requiring the invariance of the action under scale 
transformations t — > X~H, g^ — > A^ 1 ^-, one finds that 

h{X) = Cl m^X-^, f 2 (X) = c 2 m-^X^, f 3 (X) = c 3 m^X^. (107) 

The powers of X can be found from counting the scaling dimension and requiring that 
[C\ = | (note that [X] = 1); the powers of m are written to make q dimensionless. Therefore, 
instead of three arbitrary functions of X, the NLO Lagrangian for unitary Fermi gas involves 
three dimensionless phenomenological constants c±, c 2 , and c 3 (four in curved space, with 

U{x) = cM /2 x^). 

Remarkably, conformal invariance further reduces the number of independent constants 
to two (three in curved space). Performing the infinitesimal transformation (81), we find 
that Si = fCi and S4 = f £ 4 are separately invariant. Indeed, their expressions contain 
only spatial derivatives of elementary fields, and thus have conformal dimension equal to the 
scaling dimension. However, C 2 and C 3 contain time derivatives of the metric tensor and 
hence do not have well-defined conformal dimension. Thus S 2 = jC 2 and S3 = JL 3 change 
under conformal transformations, 

55 2 = 3c 2 m 1/2 J dtdx/3 [md ty /g - g ij D fl)] X 1/2 , (108) 

55 3 = -jm 3 / 2 Jdtdx^'(3X l/2 . (109) 

A closer investigation reveals that the two contributions cancel each other for a particular 
ratio of c 2 and C3. To see that, we perform integration by part in Eq. (109), and make use 
of the leading-order field equation, 

k = -t^T x + Ljt d ^ s>iD > 0)x + ^' 9 ' xd '°' (110) 

to obtain 

5S 3 = c 3 m 3 / 2 J dtdx/3 Ut^fg-X 1 ' 2 + 1-d^g g^ D 3 Q)X X I 2 + ±-Jg d.Xd.QX- 1 ' 2 

(111) 
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The contribution coming from the last term in the brackets can be further integrated by 
parts, 

J d^^g^d^d.OX- 1 ' 2 = 2 j dxy/gg^diX^djO = -2 J dxd^y/g g^d^X 1 / 2 . (112) 
One thus finds 

5S 2 + 5S 3 = (3c 2 + |) m 1 ' 2 J dt d^(3[md t ^ - g^D^X 1 ' 2 . (113) 

Therefore, conformal symmetry can be satisfied if one requires 

c 3 = -9c 2 . (114) 

(In d spatial dimensions we would have c 3 = — <i 2 c 2 .) 

Thus, we find that the NLO effective Lagrangian of unitary Fermi gas in flat space and in 
the presence of background gauge fields involves only two independent additional constants, 



£lo + L NLO = c m 3 / 2 X 5 / 2 + dy/mX-^diXdiX 

X 1 ' 2 , (115) 



+ 02 



where X — fx — D t tp - (Di<p) 2 / (2m) . 

If we retain only the trapping potential A and set A4 — O, then the effective Lagrangian 
becomes 



C = c m 3/2 X 5/2 + aVEX-^diXdiX + -^[(dtDiO) 2 - 9mV 2 A ]X 1/2 . (116) 



The two parameters c\ and c 2 are universal constants and need to be determined from 
microscopic computations. In Sec. 7 we shall see how these coefficients appear in physical 
quantities, which suggests methods to extract them. We will also show that c 2 is positive. 



7 Physical applications 

In this section, we compute some physical quantities using the NLO effective Lagrangian. 
Our aim is not to go through an exhaustive list of applications, but just to demonstrate 
the predictive powers of the EFT. Further calculations are deferred to future work. In this 
Section we use the unit system h = m = 1. 
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7.1 Phonon dispersion relation 



Consider an infinite system with no trapping potential A$ 
Lagrangian to quadratic order in <p, we find 



- ^V(Vv?) 2 + Cl (^) 2 + C2 mvV) 2 



Ai — 0. Expanding the effective 

(117) 



The phonon dispersion relation is found by solving the equation 



15Co 2 5 C 2 2, 22, 4 

— g-A*w ^ ? + Cicj g + c 2 fJ,q 



0. 



The solution is 



00 



2/i 2 

3 y 



45c 



(2ci + 3c 2 )q 4 . 



Using ji = ^€f and Eq. (83), we can write the result in the form 



uo = 



v F q 



1 



7T 



2^ ( ci + -c 2 



/up 



(118) 
(119) 

(120) 



We reproduce the linear part of the phonon dispersion with the velocity c s = a/£ /3 v p , but we 
also find the nonlinear g 3 term in the dispersion. This correction depends on the combination 
2ci + 3c 2 . In particular, if 2ci + 3c 2 < the dispersion curves up and a low-energy phonon 
can decay into two phonons. In the opposite case such a process is kinematically forbidden. 

7.2 Static density response function 

The EFT allows one to compute any correlation functions in the long-wavelength, low- 
frequency regime. Here we concentrate on the static response functions, as these functions 
are most readily computable using Monte-Carlo techniques. The static density response 
function x(k) is defined through the response of the medium to a small potential A [25]. If 
one introduces the time-ordered density- density correlation function, 



= W,x)n(0,0)), 



then 



X(0) = G nn (0,q). 

It is related to the dynamic structure function S(u, q), 

S(u,q) = -2lmG nn (uj,q), u > 0, 
through a dispersion relation 

oo 

' duj S{oj, q) 



x(q) = -2 J 



2tt 



(121) 
(122) 

(123) 
(124) 
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At zero momentum the static response function is fixed by the compressibility sum rule, 

X(0) = ~ = ~. (125) 
o/i mcj 

We shall now compute the leading correction at small q. 

From the general formalism, the density correlation function can be obtained by differ- 
entiating the generating functional W, 



G nn (t, x) 



S 2 W[A ,A t 



5A o (t,x)6A o (0,O) 



(126) 

A =Ai=0 



In EFT W can be found from diagrammatic techniques. Since to the NLO we have only tree 
graphs, W coincides with the classical action at the saddle point, 

W[A Q) A l } = S[ Vch A Q) A l ] 1 (127) 

where <p c \ is the solution to the field equation in the background A , Ai. 

To compute the static density response function one can keep Ai = and turn on only 
a static A . The leading-order field equation (110) has a trivial solution (p c \ = on this 
background. If one is interested in the NLO action one can just insert this solution into the 
NLO Lagrangian. The Lagrangian density for this solution is 



C[ Vch A ] =c (fi-A ) 5 / 2 + ^L=(d t Ao) 2 -9c 2 (V 2 A ) v ^A' , (128) 



from which we find 



x(q) 



,2 



l + 27r 2 v^( Cl -!f 2 )§f 



(129) 



At q = the compressibility sum rule is reproduced. Note that the q 2 correction is pro- 
portional to a different linear combination of C\ and c 2 compared to the phonon dispersion 
relation. 

7.3 Static transverse response function and the sign of c 2 

The transverse response function measures the "paramagnetic" current jp induced by a 
transverse static background field A: 

(j±(q)> = X T (<?)A ± (q), q • j ± = q • A ± = 0. (130) 

The "paramagnetic" current j in this formula is defined without the "diagmagnetic part" 

j = -ltytv^-V^ty), (131) 
26 



therefore if we define the transverse current- current correlator, 

i5 lk G ± (tu,q) = Jdtdxe™-*^Tj%x)j k (0,O)), i,k = x,y, (132) 
then the transverse response function is related to G± in the following way: 

X T (g)=G ± (0,g)-n. (133) 

At zero temperature, the static transverse response function approaches zero in the low 
momentum limit 

limv T (g)=0, (134) 

a fact that corresponds to the vanishing of the normal density. There is a dispersion relation 
relating x T with the dynamic transverse structure function T T (u,q): 

oo 

xr(9) = _ 2 ^r^ i (135) 



where T T is positive, so ^{q) < 0. 

The transverse response function can be computed from a formula analogous to (126), 
where A is replaced by the Aj_. Again in the presence of transverse static Ai the field 
equation yields ip = 0, and the Lagrangian density is simply 

C[A t ] = c - + c 2v ^ (-^FijFij - VdiFijA^j + 0(A 4 ± ). (136) 



We find 



X T (q) = -$c 2 Jh F q 2 . (137) 



But since x (o) < we conclude that c 2 > 0. 

7.4 Energy of unitary Fermi gas in harmonic trap 

Consider unitary Fermi gas in an asymmetric harmonic trap with the potential 

A = -{ujIx{ + u\x\ + ulxl). (138) 

To leading order, the ground state at chemical potential \i is 6 = fit, and the Lagrangian 
density, evaluated for this configuration and integrated over space, is the free energy 

-F(jm) = J c)x (co(jm - A f 2 + ^== - 9c 2 V 2 A y/^A^j . (139) 
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In Eq. (139) the integral is taken over the volume occupied by the fermion cloud, which is 
determined by A < ji. Taking the integrals one obtains 

^, , 1 u 4 /-of 9 \ uj 2 + u, 2 + Wo 9 . . 

-m = W r^ + ^\c 1 --o 1 ) J V . (140) 

where uj = (wic^^) 1 / 3 - The particle number as a function of the chemical potential is found 
by differentiating F(fi), N = —dF/d/i, and the energy as a function of particle number is 
the Legendre transform of — F: 

E(N) — F + fxN — ^(3iV) 4/3 - V2 tt 2 £ U - %c 2 \ d±A±A^ N f/\ (141) 
4 y 2 J uj 

Note that C\ and c 2 appear in the same combination as they appear in the static density 
response function (129). 

We were somewhat fortunate that all integrals in Eq. (139) converge. We would not 
be so lucky beyond the NLO. Indeed, consider as an example a 0(q 4 ) term in the effective 
Lagrangian, 

X 7/2 ■ y^ z > 

This term respects both general coordinate and conformal invariance. The contribution to 
the free energy in a trap would be 

J:: [«W ( 143 ) 



For simplicity let us assuming a spherical trap, uJi = uj. Denote the distance to the boundary 
of the fermion cloud as 5, 

x = 

then the integral in (143), in the region of small S, looks like 



6, (144) 



/ 



Clearly, the integral diverges near 5 = 0. 

This divergence can be traced to the breakdown of the derivative expansion near the 
boundary of the cloud. Indeed, near the boundary the effective chemical potential ji — 
A changes substantially over distances of order 5. When 5 is comparable to the Fermi 
wavelength, i.e., 

° r *~7^3T7«- ( 146 ) 



(/i-Ao) 1 / 2 (uj5y/ 2 n 1 / 4 ' (cj 2 /i)V6- 

the expansion parameter in the EFT becomes of order one. In other words, there is a 
boundary layer with the thickness determined by Eq. (146) where the EFT theory break 
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down. One can estimate (but not calculate) the contribution of this region to the free energy 
by placing (146) as a cutoff in (145), and find 

/x 5/3 

^boundary ~ 9/3 ' 

(147) 

The EFT, in principle, cannot give an answer with an error smaller than -^boundary This 
uncertainty translates into an uncertainty of order N 5 ^ 9 in Eq. (141). Note that this is 
smaller than the NLO correction by only a factor of N^ 1 ^ 9 . 

8 Conclusion 

In this paper we have discussed general coordinate invariance and conformal invariance 
in the nonrelativistic many-body theory. We also extract physical consequences of these 
symmetries for the unitary Fermi gas. We find that to the leading and next-to-leading 
orders in the momentum expansion we need three constants c (or £), C\ and c 2 to completely 
characterize the dynamics of the system and its response to external probes. While we have 
some numerical estimates for cq from Monte Carlo measurements of thermodynamics, we 
still lack estimates for c\ and c 2 . We propose to extract C\ and c 2 by finding the density and 
transverse response functions in the low-momentum region using Monte Carlo techniques 
and matching the results with Eqs. (7). 

The effective Lagrangian which we have derived can be used to compute many physical 
quantities. One interesting application is the computation of the kinetic coefficients of a 
unitary Fermi gas at low temperatures. The knowledge of these coefficients is important for 
computing the damping rate of oscillations of the trapped gas at low temperatures [7, 26]. 

Taking a broader perspective, we expect the general coordinate invariance to have a wide 
range of applications, since it is a common property of all nonrelativistic systems. In the 
paper we have considered only systems with one species of particles, but it is straightforward 
to generalize the general coordinate transformations (17) for systems consisting of many 
species of particles. For example, in the case of a 4 He- 3 He mixture one should introduce 
two gauge fields, one (say, A^) coupled to the 4 He current and another (B^) to the 3 He 
current. The general coordinate transformations contain two copies of Eqs. (17b) and (17c); 
in particular we have 

6A i ^--- + m 4 g ik i k , SBi = ■■■ + m 3 g ik i k , (148) 

where m 4 and m 3 are the masses of the 4 He and 3 He atoms, respectively. 

On the other hand, conformal invariance, at least in three spatial dimensions, seems to 
be specific to the unitary Fermi gas. It would be interesting to explore conformal invariance 
in lower-dimensional systems, for example for the anyon gas in two dimensions, which has 
been suggested to be a superfluid [27] and respects scale invariance [13, 28]. 
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A Christoffel symbols and Ricci tensor in the nonrelativistic limit 



Our convention for the Christoffel symbols and the Ricci tensor is 



Kx = ^ P (du9 P x + d x9pu - d p9uX ), 



R 



r> r x -i- r A v p 



The covariant derivatives of scalars and vectors are given by 



To derive Eq. (90) the following equation is used: 
For the metric of the form (27), one finds 



oo — ' v 



Y % - 
1 oo — 



m 2 c 3 



mc z 



1 0i 



AA 



2mc 2 



A ° 9a + -F s 



m 



^(V^, + V,A) + ^-, 



l 



1 °'" " 2c 9 



9k j H -^fej 

m 



(149) 
(150) 

(151) 

(152) 

(153) 
(154) 

(155) 



where and F^- are the electric and magnetic fields constructed from the field potential 

(Ao,Ai): 

Ei = d t A t - diA , Fij = diAj - djAi. (156) 

The completely spatial components of T^ x coincides with the Christoffel symbols of the 
three-dimensional space with metric gij. 
For the Ricci tensor, one has 



R oi = -7r d tdi lng + — V 



2c 



2c 



( 0w + — 



(157) 
(158) 



Again, the spatial components coincides with the Ricci tensor of the three-dimensional 
space with metric g^. In particular, in the limit c — > oo, the 4D Ricci scalar i? coincides 
with the 3D counterpart. 
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